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Abstract-This short paper studies the properties of J^-model, 
which is a human motor system model derived from equilibrium 
point hypothesis. The stability of the J^-model based on 
equilibrium and Jacobian matrix is investigated, and some 
mathematical and simulation results are presented. Especially, 
the behavior of the J^-model with tiny moment of inertia is 
discussed. The results show that the J^-model is stable and has a 
unique equilibrium point under the normal condition. But when 
the moment of inertia is tiny, the stability and equilibrium point 
will be determined by the values of the moment of inertia and 
initial points. 
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I. Introduction 

Motor disability is the one of the worst aftereffects of the 
neuro-musculo-skeletal system injury. Patients often have to 
depend on other persons or devices for even simple tasks. 
The purpose of medical management and rehabilitation 
engineering is to provide independence and mobility to a 
disabled person, which will ultimately improve his quality of 
life. However, the outstanding rehabilitation technology is 
built on the comprehensive understanding about human 
motor system. 

In this short paper, the stability of ^i-model is discussed, 
which will be useftil for deeper understanding of human 
motor system. Especially, the properties of the system under 
small moment of inertia and special descending commands 
are investigated, some mathematical and simulation results 
are presented. 

In following contents, part 2 is a concise review of the X- 
model of human motor system, which is an integrated model 
based on the equilibrium point hypothesis [1]. The equations 
of the }i-model are hsted. Part 3 is dealing with the stabilities 
of the }i-model under different situations, it suggested that the 
}i-model is stable and has a unique equilibrium point under 
certain conditions. Part 4 gives detailed investigation about 
the properties of the ^i-model when the moment of inertia is 
tiny. At last, some simulation results are presented. 

II. MODEL DESCRIPTION 

The }i-model is one of the most popular human motor 
system models, which suggests that a fixed voluntary motor 
command can be described as a set of equilibrium positions 
called an invariant characteristic. Therefore, if we change 
external torque, the limb position will shift along the 
invariant characteristic until equilibrium was once again 
established [2]. For simplicity, we selected the simple 
versions from various literatures [3,4] for simulation. 

It is assumed that the musculoskeletal system is combined 
with antagonists around a hinge joint, and the muscle 



activation may according with the changes in threshold 
angles of flexor ( /l^ ) and extensor ( A^ ), which is determined 
by the non-negative reciprocal (R) and coactivation (Q 
commands. The dynamic threshold angles ( /l^* , Aj ) depend 
on the static counterpoints and angular velocity (co = dO/ dt): 



A. =R-C,A^ = R + C,A* =A.- jUCO+ ppA* =A^- lico- p^ 



(1) 



where // is a damping factor [5]. The changes of the 
threshold angles {p^ p^ ) in equation (1) due to reciprocal 

inhibition between antagonist motor neurons adjusted by 
muscle afferents and segmental inter neurons. It is assumed 
that these inter neurons just begin their recruitment when the 
threshold is less than that for motor neurons of homonymous 
muscles, otherwise, the behavior of the inter neurons closely 
resembles that of motor neurons [3]: 

Pf =A\ -jU(0+h-9Y,p, =r[0-Aj-\-jUCO-\-hy (2) 

where r represents the strength of reciprocal inhibition, 
[xY = X if X > and [xY = if x < . 

Then, muscle activation arises when the joint angle (0) 
reaches the dynamic threshold angle, and the values are 
determined by a ftmction of the difference between them. 
Thus, Muscle activation E^ and E^ can be calculated as 
following experimental equations: 



a[e-Ar 



^aiXl-dt , 



^f-^ -\E^ = e'''^^-'' -1 (3) 

For simplicity, torques produced by the flexor and 
extensor muscles were assumed to be proportional 
to Ej and E^ , respectively: 

Tj = aEj , Te = aE^ (4) 

where a and a are constants. Equation (3) and (4) are 
consistent with experimental static torque and angle 
characteristics for constant value of X [2]. Meanwhile, 
because E^ and E^ are ftmction of velocity, they also include 

dynamic torques. The contribution of passive issues to 
elasticity and damping was ignored. The acceleration 
{£ = d^6 1 dt^) is calculated from the Newtonian equation of 
motion: 

Ie = T^-Tf (5) 

where / is the experimentally measured moment of inertia. 

III. STABILITY ANALYSIS 

Given a human motor system model, the most important 
problem is stability, because an unstable motor system is 
useless and dangerous to human. As the ^i-model is built 
based on the equilibrium point hypothesis, the determination 
of equilibrium point and stability are important. Here, we 
present new stability results of the model. The results are 
obtained based on the contraction theory presented in [6]. 
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For the ^i-model given in (1) - (5) in the previous part, it 
can be described in the state-space form [7]: 



^2 =/(^1.^2'0 



where / is given by f(x^,x^,t) = (T^-T^)/ 1, Xi and X2 are the 
state variables, corresponding to Xi=0, and X2=co, where co 
represents the angular velocity. Note that f{x^,x^,t) is a 
nonlinear piecewise function. 
Then, we have following results: 

Theorem 1: 

;; // > 0; r > 0; 7? > 0; C > 0; 

2) Xj > h. 
it has 

Xr + r/l + rh 
X+h>^ '- >X,-h 



Then, the operation regions can be combined into the 
following 5 cases: 

Xr + rX^ + rh 



7 



(1) When 

1 + r 

space equation becomes 



-<0 + jUQ) < X^-\-h, then the state- 



:[l-e 



a{xi -Af +jUX2 -rA^+rjUX2 -rh+rx^ ) 



(6) 



Thus, its equilibrium point is determined as 

Xr + rX^ + rh 



1 + r 



X,(oo)=0. 



We have the Jacobian matrix of 
1 



where 



X^h>- 



1 + r 

X^ + rXj - rh 

1 + r 



->A,-h 



a{l + r)[e 



a(x^ -ylj- +JUX2 -r/lg +rjux2 -rh+rx^ ) 



^22 = — a{iu + riu)[e 



a{xi -Xj- +//X2 -r/lg +rjiiX2 -rh+rx^ ) 



From theorem 1 , the possible operation regions of the X- 
model can be summarized to 6 cases. Consider the stability of 
system in these cases, it gives theorem 2: 

Theorem 2: For the 1-mode given in equations (l)-(5) with 
the conditions, 

1) h>0; 0<r<l; R>0; C>0; X^>h; 

// > 0; / > 0; 

2) the product axa>0. 
it has: 

rh 
for C > , there is only one equilibrium point and it is 



(2) When Xj-h<0 + juco< 



yl + rXr - rh 



/ 



1 + r 



- , then the state- 



space equation becomes 

a ^ a(A^ -jUX2 -rx^ +^^f ~^l^^2 -rh-x^ ) 



[u. 
e 
1 



(7) 



and its equilibrium point is given by 



\-r 



Xi(oo): 



X^ + rXj - rh 
1 + r 



given by 



Xr +yl 



X,(oo) = 

In this case, the Jacobian matrix has its elements 



'] 



rh 



a 



for C < , the system will not have a unique equilibrium ^22 

1-r 

point and may oscillate. 



^22 = a{iu + riu)[e 



0({Xg -fix^-rx-y +rXj- -riix^ -rh-x^) 



IV. FURTHER INVESTIGATION 



rh 

Theorem 2 shows if C < , the system may trap in 

1-r 

oscillation, but the question is whether the system should 

always oscillate under this condition? In this part, some 

investigations and simulation results are presented to provide 

a preliminary solution of this problem. 

At first, we will show the state-space equations of the system 

whenC< . IfC< , it has 



(3) When 0-\-ji(o> X^ +//, then, the state-space equation 
becomes 



1-r 



C< 



rh 



1-r 



1-r 

X^ + rXj -rh Xj -\- rX^ + rh 



In this case, the equilibrium point is 

(x^(oo) = X^ 
[X,(oo) = 

and the Jacobian matrix has the form 
1 



(8) 



1 + r 



1 + r 
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(4) When 6 + juco <Aj-h, then the state-space equation 
becomes 



i^=^ = ^[e-(^^-/'^^-i)-l] 



In this case, its equihbrium point is determined as 

X,(oo) = 

and the Jacobian matrix has the form 
1 



-a[e 



a(A^-jUX2-x^) 



-aju[e 



a(A^-jUX2-x^) 



A + rAr -rh Xr-V rX^ + rh 
(5) When^ ^ <e + jU(o<^—^ , then, the 



1 + r 
state-space equation becomes 

and the Jacobian matrix has the form 

"0 r 



1 + r 



(10) 





Then we get the operation regions in 6-cd space as shown 
in Fig. 1. Note that the system states with the initial points 
given in items (l)-(4) may enter the adjacent regions as their 
equilibrium points are in the adjacent regions or on borders. 
But the middle region in item (5) is an unstable region. Thus, 
if the system converges monotonically to the equilibrium 
points, then it may have two equilibrium points as in items 
(1) and (2). If the system converges to the equilibrium points 
with oscillatory behavior, then it may enter the middle region. 
In this case, it will cross over the middle region and enter the 
next region opposite to the previous region. Therefore, the 
system may oscillate. 




{Xfir\^ky{\+?^<e^,oi<\+h 



Fig.l: Operation regions when C < 



rh 



Therefore, the problem at the start of this part becomes: 
how can we judge whether the system will converge 
monotonically to the equilibrium point or not. Because the 
system states with initial points given in items (3) and (4) 
must enter the adjacent regions, we will only consider the 



regions defined in items (1) and (2), that is, the regions 
border upon middle region. 

As mentioned before, the general format of the ^i-model's 

1 

, therefore, the eigenvalues are: 



(9) Jacobian matrix is 



4 =-K2± 74^ + ^22'] 

According to the knowledge of Jacobian matrix, if the 
eigenvalues have the imaginary part, the system will trap in 
oscillation. Therefore, it gives 

4^21 + ^2/ < 

if the system converge to the equilibrium point with 
oscillatory behavior. 

Apparently, in items 2 and 3, it has a22^l^ci2\ and <22i<0. 
Then we have: 

4^21 + ^22^ = (4 + ji^a^^a^^ < ^ 4 + ji^a^^ > ^ jd^a^^ > -4 
Combined with expressions of items (1) and (2), the 

above inequation becomes: 
item (1): 

2 

a(x^-Af+jUX2-rA^+rjUX2-rh+rx^)-. . 



2 ^J^ /1 xr ' 

// ^21 = — ^a(l-\-r)[e 



item (2): 



^^a,,=-^^a{\ + r)[e'"''-''^- 



+rAr-rf^X2 -rh-x^ ) 



]>-4 



^^a(\ + ^)r^^(^^-^^2-^^i+^^/-^/^X2-r/z-Xi) -. ^ ^ 



Therefore, if the moment of inertia is large enough, these 
inequations will be proved, that is, the system states with 
initial points given in items (l)-(4) will oscillate around 
middle region; otherwise, the system will converge 
monotonically to the equilibrium points on the border of 
middle region {Xf and Xe )• This result suggests that if the 
moment of inertia of human limb is too small, the limb 
movement will stop stably at an equilibrium position, which 
is determined by the initial states. 

IV. SIMULATION RESULTS 
For verifying the above mathematical conclusions, 
simulation is carried out. 

Model is used for simulation with following parameters 
[3]: 

7^=60% C=r,r=0.2, 

7=0.115 kg.m', a=1.2N.m, 6^^=0.05 deg\ /?=10°,//=0.58 s. 
Then, the threshold angles of fiexor {Xj) and extensor {X^) 
can be calculated as: 
X=R-C=59\X=R + C=6\\ 
rh 0.2*10 



and 



:2.5°>C 



\-r 0.8 

At first, the performance of this model under different 
moments of inertia is investigated. Fig. 2 shows the system 
dynamic responses with different initial values and moments 
of inertia. The blue curve and red curve start fi-om same 
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initial point with the different moments of inertia. The result 
shows that the blue curve oscillates around the middle region, 
but red curve ceases as a horizontal line, which is also the 
upper border of middle region. It suggests that the movement 
of the limb with small moment of inertia will stop at a certain 
position, even under the impotent descending commands. 
Note the pink curve and red curve start from different initial 
points with the same moment of inertia. These two curves 
converge to different horizontal lines. Because the system 
enter the middle region with different velocity, and the 
acceleration of the middle region is zero, that is, the system 
will keep its own velocity before entering the next region, 
and will stop at the border of next region by the enormous 
deceleration. 




Time (s) 



Fig. 2. The dynamic response of the system with different initial points (9, co) 
and moment of inertia (7) 



1=0.00115 kg*m2 

R=60 degrees 

C=1 degrees 

|j=0.58 s 

r=0.2 

h=10 degrees 

Time: 0- 1000 s 




Joint Angle (degree) 



Fig. 3. The dynamic response of system with different initial points (0, co) 
and the tiny moment of inertia. 

Then we face a new problem: when the system will enter 
the middle region and converge to the border of opposite 
region, and when the system will monotonically converge to 
the equlibrium point at the border without entering the middle 
region. Fig 3 shows the dynamic response of the system with 
different initial points and the tiny moment of inertia. It 
suggests that the system states near the two candidate 
equilibrium points determine the fmal equilibrium point. If 
the system tends to the equilibrium point of current region 
before entering the middle region, and is decelerated when it 



is close to this equilibrium point, the system will not span the 
middle region. Otherwise, the system will converge to the 
equilibrium point in opposite region. 

V. Conclusion 

Human motor performance appears to be remarkably 
flexible and dexterous; this requires the model of human 
motor system must be stable and easy to control. Therefore, 
the discussion for properties of the model under different 
conditions is very necessary and significant. In this short 
report, an important human motor system model - }i-model is 
described and investigated. The results about stability of this 
model are presented, which suggest that the }i-model will be 
stable and have a unique equilibrium point under certain 
conditions. Especially, when the moment of inertia of limb is 
tiny, the stability and fmal equilibrium point will be 
determined by the values of moments of inertia and initial 
points. These results should be served as a guideline on the 
forecast for the behavior of some special motor systems, and 
may be usefril for studying the stability of other nonlinear 
piecewise frinctions. 
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